Limits on the WIMP-nucleon interactions with CsI(Tl) crystal detectors 
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The Korea Invisible Mass Search (KIMS) experiment presents new hmits on the WIMP-nucleon 
cross section using data from an exposure of 3409 kg-d taken with low-background CsI(Tl) crystals 
at Yangyang Underground Laboratory. The most stringent limit on the spin-dependent interaction 
for a pure proton case is obtained. The DAMA signal region for both spin-independent and spin- 
dependent interactions for the WIMP masses greater than 20 GeY/c^ is excluded by the single 
experiment with crystal scintillators. 

PACS numbers: 95.35.-|-d, 14.80.Ly 



The existence of dark matter has been widely sup- 
ported by many astronomical observations on vari- 
ous scales [1][2] 3]. Weakly interacting massive parti- 
cles (WIMPs) are a good candidate for dark matter well 
motivated by cosmology and supersymmetric models ^3] . 
The Korea Invisible Mass Search (KIMS) experiment has 
developed low-background CsI(Tl) crystals to detect the 
signals from the elastic scattering of WIMP off the nu- 
cleus [ilifT!]. Both Cs and ^^'^I are sensitive to the 
spin-independent (SI) and spin-dependent (SD) interac- 
tions of WIMPs. Recently, the role of Csl in the direct 
search for SD WIMP for pure proton coupling has been 
pointed out It is worth noting that ^-^^I is the dom- 
inant target for the SI interactions in the DAMA exper- 
iment. The pulse shape discrimination (PSD) technique 
allows us to statistically separate nuclear recoil (NR) sig- 
nals of WIMP interactions from the electron recoil (ER) 
signals due to the gamma ray background 0] [lol | . 

The KIMS experiment is located at the Yangyang Un- 
deground Laboratory (Y2L) at a depth of 700 m under an 
earth overburden. Details of the KIMS experiment and 
the first limit with 237 kg-d exposure data can be found 
in the previous publication [ll|. Four low-background 
CsI(Tl) crystals are installed in the Y2L and operated 
at a temperature of T = 0°C. Throughout the exposure 
period, the temperature of the detector was kept sta- 
ble to within ±0.1°C. Green-enhanced photomultiplier 
tubes (PMTs) are mounted at both ends of each crystal. 
The signals from the PMTs are amplified and recorded 
by a 500 MHz FADC. Each event is recorded for a pe- 
riod of 32 ^s. Both PMTs on each crystal must have at 
least two photoelectrons within a 2 /is window to form an 
event trigger. We obtained 3409 kg-d WIMP search data 



TABLE I: Crystals used in this analysis and amount of data 
for each crystal 



Crystal 


mass (kg) 


data (kg-days) 


S0501A 


8.7 


1147 


S0501B 


8.7 


1030 


B0510A 


8.7 


616 


B0510B 


8.7 


616 


Total 


34.8 


3409 



with four crystals, as shown in Tablc[Tl The energy is cali- 
brated using 59.5 keV gamma rays from an Am source. 
For calibration of the mean time, a variable used for the 
PSD, NR events are obtained with small crystals ( 3 cm x 
3 cm X 3 cm ) using an Am-Be neutron source. Compton 
scattering events taken with the WIMP search crystals 
using the ^^^f^g source are used to determine the mean 
time distribution of the gamma background. Compton 
scattering events are also taken with the small crystals 
to verify that the mean time ditributions for both the test 
crystals and the WIMP search crystals are the same. In 
order to understand the nature of the PMT background, 
a dominant background at low energies, acrylic boxes are 
mounted on the same PMTs used for the crystals. The 
data obtained using this setup is used to develop the cuts 
for the rejection of PMT background. 

Since the decay time of the scintillation light in the 
CsI(Tl) crystal is rather long, photoelectrons are well 
separated at low energies and thereby enabling recon- 
struction of each photoelectron. The time distribution of 
photoelectrons in an event is fitted to a double exponen- 
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FIG. 1: (color online). MT distribution of NR events (open 
squares), ER events (open circles) and WIMP search 
data (filled triangles) of S0501A crystal in the 5-6 keV range. 
Fitted PDF functions are overlayed. x^/DOF =0.8 and 1.3 
with DOF=3S and 35 for NR and ER events respectively. 
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FIG. 2: (color online). Extracted NR event rates of the 
S0501A (open circles), S0501B (filled circles), B0510A (filled 
squares), and B05010B (filled triangles) crystals and only sta- 
tistical errors (la) are shown. The points are shifted with 
respect to each other on the 2;-axis to avoid overlapping. 



tial function given by 

/(t) = _ exp <^ 



R 

— exp 



where r/ and Ts are decay time constants of fast and 
slow components, respectively, R is ratio between two 
components, and to is the time of the first photoelectron 
in the event. The mean time (MT) of each event is then 
calculated using these quantities as 



t ■ f{t)dt/ / fit)dt 



MT : 



With this method, an improvement in PSD is achieved 
over the previous analysis where we used a simple math- 
ematical mean In order to reject the PMT back- 



ground, we applied cuts to the fit variable, t/ . The ratio 
between the maximum log likelihood value of the dou- 
ble exponential fit and that of the single exponential fit 
is also used to reject the PMT background, since PMT 
background events tend to be shaped as single exponen- 
tial decay. To reject the background that originates from 
the radioactivity of the PMT, the asymmetry between 
the signals from two PMTs is applied. Finally events in 
which signals are recorded in more than one crystal are 
rejected. The event selection efficiency was estimated 
by applying the same analysis cuts to the neutron and 
gamma calibration samples. The efficiency depends on 
the measured energy and ranges from 30% at 3 keV to 
60% above 5 keV. 

The estimation of the NR event rate is performed in 
each 1 keV bin from 3 to 11 keV for each crystal. The MT 
distributions of NR events and ER events are compared 
with the WIMP search data in Fig. [1] for the 5-6 keV 
energy range. The probability density functions (PDF) 
for the ER and NR events are obtained by fitting these 
distributions. An unbinned maximum likelihood fit is 



performed with the log{MT) distribution of the WIMP 
search data using the likelihood function, 

n 
k=l 

where the index i denotes the i-th energy bin; n — 
NNR,i + NER,i is the total number of events; N^R^i and 
Ner^i are the numbers of NR and ER events, respec- 
tively; PDFMR,i and PDFer,, are PDFs of NR and ER 
events, respectively; and = log{MT) for each event. 
The NR event rates obtained for each bin and for each 
crystal after efficiency correction are shown in Fig. [2l 
The extracted NR event rates are consistent with a null 
observation of the WIMP signal. 

In order to obtain the expected measured energy spec- 
trum of a WIMP signal including instrumental effects, 
a Monte Carlo (MC) simulation with GEANT4 [13] is 
used. A recoil energy spectrum is generated for each 
WIMP mass with the differential cross section, form fac- 
tor, and quenching factor, as described in Ref. [ll]. The 
spin-dependent form factor for ^^'^Cs calculated by Toiva- 
nen [l4| is used, while for ^^^I, Ressell and Dean's cal- 
culation [l^ is used. The photons generated with the 
fitted decay function described above are propagated to 
the PMT and digitized in the same manner as in the ex- 
periment. Subsequently, the photoelectrons within given 
time windows are counted to check the trigger condition 
and to calculate energy. In this manner, the trigger ef- 
ficieny and energy resolution is accounted for in the ex- 
pected energy spectrum. The trigger efficiency is found 
to be higher than 99% above 3 keV. The simulation is 
verified with the energy spectrum obtained using 59.5 
keV gamma rays from ^'^^Am. The peak position and 
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TABLE 11: Spin expectation values for "^Cs and ^^''l 



Isotope 


J 


< Sp> 


< S„ > 


Reference 




7/2 


-0.370 


0.003 


fl6] 


127j 


5/2 


0.309 


0.075 


fl5] 



width of the distribution are very weU reproduced for 
each crystal as described in Ref [11,]. 

The total WIMP rate, R, for each WIMP mass is ob- 
tained by fitting the measured energy spectrum to the 
simulated one. The 90% confidence level (CL) limit on 
R is calculated by the Feldman-Cousins's approach in 
the case of Gaussian with a boundary at the origin [l3| 
and then converted to the WIMP-nucleus cross section, 
aw- A ■ Subsequently, the limits on WIMP-nucleon cross 
section is obtained from Ref. l3| 18 1 as follows: 



aw- 



where fin,A are the reduced masses of the WIMP-nucleon 
and WIMP-target nucleus of mass number A. Ca/Cu — 
for SI interactions and Ca/C„ — 4/3{ap < Sp > 
+a„ < Sn >}^iJ + ^)/J for SD interactions. Here Op, 
a„ are WIMP-proton and WIMP-neutron SD couplings 
respectively. The spin expectation values used for this 
analysis are shown in Table |TT1 Following the "model- 
independent" framework 18], we report the allowed re- 
gion in two cases for SD interaction: one for a„ = 0, 
and the other for Up = 0. We express the WIMP-nucleon 
cross section as follows: 



1 



a 



W—n,p 



/4 _ 

aW~A 



J 



H\ 4 (J+ 1) < Sn,p 



where we indicate pure proton (p, a„ = 0) and pure 
neutron (n, Op = 0) coupling for SD interaction. We also 
present the allowed region in the Up — a„ plane with the 
following relation [18]: 
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where Gp is the Fermi coupling constant. 

The uncertainty in the MT distribution results in the 
uncertainty of the NR event rate. The limited statistics 
of the calibration data and different crystals used for the 
neutron calibration and WIMP search data are the ma- 
jor sources of this uncertainty. The former is investigated 
by varying the fitted parameters in PDF function within 
errors. The lattter is estimated by changing the mean of 
MT by the difference between the crystals. The system- 
atic uncertainties from these two souces are combined in 
quadrature resulting in 20-30% of statistical uncertain- 
ties depending on the energy bins. In addition, there 
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FIG. 3: (color online). Exclusion plot for the SD interaction 
in the case of pure proton coupling (a„ = 0) at the 90% 
confidence level 



lO-" 



10= 



a. 
c 

o 

0) 

<n 

I 10 

o 

i_ 

o 
a 

CO 

3 

0) 

Z 

1 10-' 
i 

10 = 




lO'' 10^ 
WIMP Mass (GeV) 



10* 



FIG. 4: (color online). Exclusion plot for for the SD interac- 
tion in the case of pure neutron coupling (a^ = 0) at the 90% 
confidence level 

are uncertainties in the MC estimation of the expected 
event rates due to the uncertainties in the quenching fac- 
tors and the difference of energy resolution between the 
MC simulation and the data. The systematic error from 
the MC simulation is estimated to be 13.3% of the limits. 
These systematic errors are combined with the statistical 
error in quadrature in the presented results. 

The limits on the SD interactions are shown in Fig. [3] 
and [4] in the cases of pure proton coupling and pure neu- 
tron coupling, respectively. We also show the results ob- 
tained from CDMS NAIAD 20], SIMPLE [U, and 
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FIG. 5: (color online). Allowed region (90% confidence level) 
in Bp — a„ plane by KIMS data (inside the solid line contour) 
for 50 GeV WIMP mass. Results of CDMS Igj (dotted line) 
and NAIAD [20l| f dot-dashed line) are also shown. 



masses higher than 20 GeV/c^ by the single experiment. 
The most stringent limit on the SD interaction in the 
case of purely WIMP-proton coupling is obtained. 
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FIG. 6: (color online). Exclusion plot for the SI interactions 
at the 90% confidence level. 

PICASSO dil. The DAMA signal region is taken from 
Ref ^23'!. Our limit provides the lowest bound on the 
SD interactions in the case of pure proton coupling for a 
WIMP mass greater than 30 GeV/c^. The allowed region 
in the Op — a„ plane for the WIMP mass of 50 GeV/c^ is 
also shown in Fig. [5] together with the limits from CDMS 
and NAIAD. The limit for the SI interactions is shown 
together with the results of CDMS EDEL- 
251, CRESST 26], ZEPLIN I [l^, and the 3cr 



in Fig 
WEISS 

signal region of DAMA (1-4) [28|. Although there are 
several experiments that reject the DAMA signal region, 
this is the first time that it is ruled out by a crystal de- 
tector containing ^^^I, which is the dominant nucleus for 
the SI interactions in the Nal(Tl) crystal. 

In summary, we report new limits on the WIMP- 
nucleon cross section with CsI(Tl) crystal detectors using 
3409 kg-d exposure data. The DAMA signal regions for 
both SI and SD interactions are excluded for the WIMP 
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